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Abstract 

We study instanton contribution to the partition function of the one matrix model 
in the k-th multicritical region which corresponds to the (2, 2k — 1) minimal model 
coupled to Liouville theory. The instantons in the one matrix model are given by local 
extrema of the effective potential for a matrix eigenvalue and identified with the ZZ 
branes in Liouville theory. We show that the 2-instanton contribution in the partition 
function is universal as well as the 1-instanton contribution and that the connected 
part of the 2-instanton contribution reproduces the annulus amplitudes between the 
ZZ branes in Liouville theory. Our result serves as another nontrivial check on the 
correspondence between the instantons in the one matrix model and the ZZ branes in 
Liouville theory, and also suggests that the expansion of the partition function in terms 
of the instanton numbers are universal and gives systematically ZZ brane amplitudes 
in Liouville theory. 
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1 Introduction 



The nonperturbative study of noncritical strings in terms of the matrix models 12 0] led 
to an important suggestion that string theory in general possesses the nonperturbative 
effect that behaves as ~ e~9^. It was pointed out in 5J and is now widely believed that 
this effect is attributed to D-branes. Indeed, the discovery of D-branes in Liouville theory 
which are called the FZZT branes jSllZllH] and the ZZ branes [0] triggered recent progress in 
noncritical strings and the matrix models, in which the origin of the nonperturbative effect 
e~9i given by the string equations of the matrix models jH ^21 ^] was identified with the 
ZZ branes P]-[T?)j. 

The authors of Ref. J7j studied intensively the D-branes in a series of noncritical strings, 
the (p, q) minimal conformal field theory coupled to two-dimensional quantum gravity (Li- 
ouville theory).^ They call this the {p,q) minimal string theory. They showed that there 
exist '^P"-^)^^^"-'-^ independent ZZ branes labeled by (m, n), where qm — pn > 0, and those ZZ 
branes correspond to the singularities of an auxiliary Riemann surface that are formed by 
the analytically continued boundary cosmological constant and the derivative of the FZZT 
disk amplitude with respect to the boundary cosmological constant. In particular, when 
{p, q) = (2, 2k — 1), the minimal string theory is realized by the A;-th multicritical region of 
the one matrix model ^Hl and the ZZ branes are identified with the local extrema of the ef- 
fective potential for a matrix eigenvalue ^71^]. These extrema can be called the instantons. 
The annulus amplitudes between the D-branes in the minimal string theory were evaluated 
in |2n| to examine the deformations by the D-branes. 

The authors of Ref. explored in detail the nonperturbative effect stemming from the 
ZZ brane in c = noncritical string theory (the (2, 3) minimal string theory) from the 
viewpoints of the one matrix model as well as of the loop equations (the string field theory). 
One of the important results of Ref. [TU] is that the ratio of the 1-instanton sector in the 
partition function of the matrix model to the 0-instanton sector is universal, namely, it does 
not depend on the detailed structure in the potential of the matrix model. Actually, they 
confirmed this up to next to leading order in the expansion by the explicit calculation. 
This ratio is interpreted as the chemical potential of the instanton. The results of and 
[in] tell us that the leading order of this ratio is equal to e^^^, where Zzz is the ZZ brane 

^Here, p and q are relatively prime integers and p < q. 
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disk amplitude and Zzz ~ ^1 9s- The extension of this analysis to the supersymmetric case 
(c = type OB string theory) was reported very recently j21j . 

In general, the interaction between D-branes is given by the annulus amplitude between 
the D-branes. If the identification of the instantons in the matrix model with the ZZ branes 
in Liouville theory is valid, we should be able to obtain the annulus amplitudes between 
the ZZ branes by considering the interaction between the instantons. This means that we 
go beyond the dilute gas approximation. For this purpose, we need to consider the general 
(2, 2A; — 1) case. The reason is as follows. The interaction between the identical instantons 
diverges due to the Vandermonde determinant for the matrix eigenvalues. Correspondingly, 
the annulus amplitudes between the identical ZZ branes also diverges. On the other hand, 
the annulus amplitudes between the different ZZ branes is finite. The (2, 2k — 1) theory 
possesses k — \ independent instantons in the one matrix model and the k — \ independent 
ZZ branes in Liouville theory. So, the (2, 3) case is not sufficient for our purpose. 

First, we generalize the calculation of the ratio of the 1-instanton sector to the 0-instanton 
sector in JHl to the (2, 2k — 1) case. We show that it is also universal for generic k and the 
leading order is given by the ZZ brane disk amplitude. Next, we consider the 2-instanton 
sector in the partition function of the matrix model turning on the interaction between the 
instantons. We show that the ratio of the 2-instanton sector to the 0-instanton sector is 
universal and its connected part indeed reproduces the annulus amplitudes between the ZZ 
branes in Liouville theory. Our result serves as another nontrivial check on the correspon- 
dence between the instantons in the one matrix model and the ZZ branes in Liouville theory, 
and also suggests that the expansion of the partition function in terms of the instanton 
numbers are universal and gives systematically the ZZ brane amplitudes in Liouville theory. 

This paper is organized as follows. Sections 2 and 3 are devoted to the reviews of the 
ZZ branes in the (2, 2k — 1) minimal string theory and of the k-th multicritical region of the 
one matrix model, respectively. In section 4, we see the behavior of the effective potential 
of a matrix eigenvalue and define the expansion of the partition function in the instanton 
numbers. Sections 5 and 6 are the main part of this paper. We perform the above mentioned 
calculations of the ratio of the 1-instanton sector to the 0-instanton sector in section 5 and of 
the ratio of the 2-instanton sector to the 0-instanton sector in section 6. Section 7 is devoted 
to summary and discussion. In appendix A, we give the detailed calculation for the (2, 5) 
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case. In appendix B, we gather some formulae used in the main text. 



2 ZZ branes in Liouville theory 

In this section, we describe a part of the results in Refs. pT| ED] which is relevant for our 
purpose. Many of the features of the {p, q) minimal string theory are provided by an auxiliary 
Riemann surface M.p,q-, which is described by the algebraic equation 

F(e,r/)=r,(O-T,(r;)=0, (2.1) 

where Tp(cos^) = cospd is the Chebyshev polynomial of the first kind. ^ is the ratio of the 
boundary cosmological constant ( to the square root of the bulk cosmological constant fi 
while 1] is proportional to the derivative of the FZZT disk amplitude Zpzzr with respect to 
C: 

C 

{=— , ?7~/i ^pdt^ZpzzT- (2.2) 
It is convenient to introduce auxiliary parameters a and z = cosh in terms of which 

^ = cosh7r,^^cr = Tp{z), 

1] = cosh TT A -a = T„(z). (2.3) 

V p 

Note that z covers the surface exactly once. The ZZ branes correspond to the singularities 
of A^p,g given by F = d^^F = dy^F = 0, which correspond to two different values of z denoted 
by z"^. 

In what follows, we restrict ourselves to the case in which {p,q) = {2,2k — 1). In this 
case fl2.1|) is given by 



= T..-i(0 + 1 = ^^^^^^+^^^^^^^^^\ (2.4) 



We define ?7a:(0 by 



v/2,»(f) = (2.5) 



Then, F = d^F = drjF = are equivalent to 

Vk = 0, %fc(e)2 = 0. (2.6) 



These equations are solved as 



which corresponds to 



27rn 

C„ = -cos^^— y, n = 1,2, (2.7) 



, ivrn 
^„=-sin^^. (2i 



characterizes the (l,n) ZZ brane. r^^ is expressed in terms of ^^'s as 

2%(0^ = 2''~'{^ - - 6)' ■ ■ • - + 1). (2.9) 

It is convenient to introduce an integral of r/k- 



y ^^^^^ ^ = (2fc+l)v/eTT - (2^-3)v/m- ^'-'^^ 



This is proportional to the FZZT disk amplitude: ~ A* ^ ^ZpzzT- Later, we will use 

the following quantities which are proportional to the ZZ brane disk amplitudes. 

..te.)^(-l)-v^(^ + ^)si„J!^. (2.11) 

The annulus amplitudes between the ZZ branes were calculated in The result for that 
between the (l,Ti) and (l,Ti') ZZ branes is 

3 The /c-th multicritical region of the one matrix model 

We are concerned with the one matrix model with a generic potential. 

^ oo 

V{x) = -x'-Y. (3-1) 

m=3 

where is an x Hermitian matrix. By diagonalizing 0, this integral is reduced to 

N 



Z= I nrfA,A;v(Ai,---,A^)V^S£i 
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(3.2) 



where Ai, ■ ■ ■ , Aat are eigenvalues of and AAr(Ai, ■ ■ ■ , Xn) is the Vandermonde determinant 
in terms of Ai, ■ ■ ■ , Ajv- It is a standard technique to introduce the orthogonal polynomials 
P„(x), which satisfy 

J dx e-^^(^)p^(x)P„(x) = hJmn, (3.3) 

where P„(x) is a polynomial of degree n and is normalized so that the coefficient of equals 
one. It is easy to see that the following recursion relation holds. 

), (3.4) 

Tn = T^. (3.5) 

The partition function is expressed in terms of r„: 



N~l 



Z = Nlhoh, ■ ■ ■ hN-i = N\h^ Yl r^-'^. (3.6) 



n=l 



Then, the relevant part of the free energy F = log Z takes the form 



7V-1 



F = 5^(iV-n)logr„. (3.7) 



n=l 



The Schwinger-Dyson equations for the orthogonal polynomials, 

/ ax 



= Jdx A(e-^^Wp„(x)P„(a;)), (3.8) 



give recursion relations for r„ and from which one can determine and s„ as functions 
of g^s. 

We need the k-th multicritical region of the one matrix model to obtain the (2, 2k — 1) 
minimal string theory. The k-th multicritical region is realized by fine-tuning k—1 parameters 
among g^s and taking — oo limit. We introduce a continuous variable cr = in order 
to examine the critical behavior of the model. The critical point corresponds to a = 1. As 
is explained in appendix C of Ref.jinj, if the functions r(cr) and s{a) are defined by 

rn = r{a), 

Sn = ^(^ + ^)' (3-9) 



r(cr) and s{(t) are 0{N^) quantities and the corrections start with 0{^). Moreover, the 
(9(A^°) parts of r(cr) and s{a) behave at the critical point hke 



Or 
da 
ds_ 
da 



(3.10) 



where Vc is the critical value of r^. Now we are ready to write down the scaling limit of 
the one matrix model which gives rise to a perturbation around the k-th critical point and 
corresponds to the (2, 2k — 1) minimal string theory 



9r, 



■,k-l, 



r(cr) = Tc 



-aeuir) 
2 ^ ' 



s{a) = Sc- -ay/r'ceu{T) 

a = 1 — e^ur, 
1 

N 



k4-l 



(3.11) 



where e is a cutoff so that e — corresponds to the continuum limit, grmci fc and Sc are 
critical values of grm, r and s, respectively, which are dependent on the detailed structure 
in the potential of the matrix model, /i is the bulk cosmological constant which is identified 
with that in the Liouville theory, a, (3rni, ^ and k are certain constants. In ()3.11|) . we restrict 
ourselves to the leading order of the 1/A^ expansion and have taken ()3.1()|1 into account, a, 
/3m, and v in ()3.11|1 are adjusted in such a way that u{a) obeys a string equation [H 



(3.12) 



where 



other tj = 0. 



P 



0,1,- 



TT 



2^— 2p 



v/8 ik-2p)\p\T{p-k + l)' 



(3.13) 



This represents the above mentioned perturbation around the k-th multicritical point. When 
r = 0, the string equation ()3.12|) allows a solution 



u{0) = y/JI. 



(3.14) 



The universal part of the sphere contribution to the free energy ()3.7|) is expressed by u{t) 

as 

pisphere) ^ j da (1 - a) logr{a) = i^"^ ag^^ J drruir). (3.15) 

In the following sections, we use the resolvent, which is defined in the large limit (the 
leading order of the expansion) by 

i?(x) = /i-tr f^-^y (3.16) 



R{x) is related to the eigenvalue density p(x) as 



p(x) = --lmR{x + iO). (3.17) 
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By solving the loop equation, the form of R{x) is determined as 

R{x) = ^V'{x) + W{x), 

W{x) = iv^'(a;)2+p(x), (3.18) 

where p{x) is a polynomial of degree mo — 2 when 

gmo 7^ 0, gm = for m> itlq. (3.19) 

p{x) is determined by the structure of the cut, namely the location where the eigenvalues 
are distributed, and the condition that R{x) ~ ^ when |x| — > oo. We are interested in the 
one-cut solution, in which p[x) is nonzero only for the period [6, a]. Then W{x) takes the 
form 

W{x) = -K{x){x — Xi){x — X2) ■ ■ ■ (x — Xk-i)^/ {x — a)(x — b), (3.20) 

where K{x) is a polynomial of degree mo — k — 1. In the scaling limit ()3.1H) . a and x„ behave 
like 

a = X*(l - XaA//i^), 

Xn = x^{l + Xn^/Jj^), n = 1,- 1, (3.21) 

where Xa and Xn are given constants. On the other hand, one can regard b and K{x) as 
some constants in the scaling limit. If x is scaled in the scaling limit ()3.11|) as 

x = x,{l + eC), (3.22) 



W{x) starts with a term proportional to e*^ 2 which is the universal part of the resolvent. We 

1 

will see in section 5 that Xn = (.nXa and Xa = ax'^Vc for generic V{x), Xa = jCi for even V{x). 
Hence, if ( is tuned as C = XaC? the universal part of the resolvent becomes proportional to the 
derivative of the FZZT disk amplitude with respect to the boundary cosmological constant, 
namely ri{^). This is anticipated because the resolvent is interpreted as the expectation value 
of a marked macroscopic loop in the matrix model and the macroscopic loop is nothing but 
the geometrical meaning of the FZZT brane. In appendix A, we illustrate the calculations in 
this section with the case in which k = 3 and the potential is even. We see that Xn = inXa 
and Xa = \(^ actually hold. 

4 The effective potential for an eigenvalue and instan- 
tons 

We consider the situation in which a single eigenvalue, say Aiv, is separated from the others. 
The partition function ()3.2p is expressed as 



where we set Aat = x. By using an (A^ — 1) x (A^ — 1) Hermitian matrix (pN-i, this is rewritten 

as 

Z = Z^_iy"cix(det(x-0^_i)2)jv.ie-^^(^), (4.2) 



where 



Zn_i = / e 



~NtTV(<t,N-i) 



{0)m-i = / ti0;v-i O e-^t^^(^--) (4.3) 



The effective potential for x is defined by 



Veffix) = V{x) - ^log(det(x - (j)N-iy)N-i (4.4) 



in such a way that 

Z = ZM-ij dx e-^^=//(^). 



(4.5) 



At leading order of the expansion, the following calculation is justified. 
(det(x 



^n-iY)n-i 



(det(x-0)2) 
exp [2Re(trlog(x — 0))] 



exp 



/X 
dx R{x') 



(4.6) 



Therefore, using ()3.18|1 . we find that the leading order of Ve//(x) is 



/X 
dx'W{x'). 



(4.7) 



Or equivalently, V 



(0), 



-2ReW{x). Then, we see from flH.2()|l that v}fl{x) is constant in 



the cut and that V^jj-{x) takes local extrema at x = xi, ■ ■ ■ , Xk-i- We ignore extrema coming 
from K{x) in ()3.20|) . since they do not contribute in the scaling limit. In Fig. 1, we draw 
the shape of V}jj{x) roughly in the k = 5 case. 







\ X2 / 








b 


a xi X3 X4 X 



Figure 1: Effective potential for an eigenvalue in the k = 5 case 

These extrema can be considered as the instantons in the one matrix model and will be 
identified with the ZZ branes in Liouville theory. We label by {gi, ■ ■ ■ , qt-i} the configuration 
in which qn {n = 1, ■ ■ ■ , k — 1) eigenvalues among are located around x = Xn and the other 
N — q eigenvalues are located in the cut, where q = qi + ■ ■ ■ + g^-i- Namely, g„ is the 
instanton number of the n-th instanton. We denote by dx the perturbative expansion 
around the 'classical solution' x = Xn, which yields a expansion. The leading and 
subleading contributions to this expansion are nothing but the saddle point integral over 
X around x = Xn- We expand the partition function in terms of the instanton numbers as 
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follows. 

oo 

2i — ^ ^ ^{gi,—,gk-i} 

9l,---,9fc-l=0 

oo 

_ ^(O-inst.) J^{qi,-,qk-i} ^ 

qi,---,qk-i=o 

where 



gi! ■ ■ ■ (A^ - q) ! 1^ .^j-^ y ib<A.<a 

fc-1 gn N-q 



X 



n n n - ^^)' ^^-"^(^i' • • • > ^^-'^)' 



Vn=l i™=l i=l 



r ^ 

/l{9i,--->?fc-i} — 

_^(0-mst.) ' 

^{°'-'°> = 1. (4.9) 

For example, 

_2-{0,-,0,9n=l,0,-,0,q„,=l,0,-,0} 



iV-2 



N{N -1) f dx f dy f Yl 

J Xn J^n' Jb<Xi<a 

/N-2 



We make the following abbreviations for the quantities which we are concerned with in 
subsequent sections. 

v{0,-,0,g„=l,0,-,0} 
^{0,-,0,q„=l,0,-,0} ^ ^ A{n) 

7{0,-,0,g„=2,0,-,0} 

^{0,-,0,g„=2,0,-,0} ^ ^ ^(n,n) 

ZiO-inst.) ' 

7{0,-,0,(j„=l,0,-,0,g„/=l,0,-,0} 

^{0,-,0,9n=lA-,0,q„,=l,0,-,0} ^ ^ ^(nX) ("4 J^]^) 

Z{0—inst.) ' V • / 
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A^^'' is the ratio of the 1-instanton sector to the 0-instanton sector while ^("'") and A^"^'""'^ 
are the ratios of the 2-instanton sectors to the 0-instanton sector. We make ^("'"') represent 
both the cases, n — n' and n ^ n' . The 1/N expansion on which our calculations in the 
following sections are based is not the expansion in term of l/N"^ but the one in terms of 1/N 
due to the instanton effects. It will turn out that log^^^^ and log^^"'"') start with 0{N). 
We will evaluate 0{N), 0{\ogN) and 0{N^) terms, which become 0{\oggs) and 

0{g1) terms in the continuum limit, respectively. 

5 1-instanton sectors and the ZZ brane disk ampli- 
tudes 

In this section, we calculate A^'^\ log^^"^ will turn out to start with 0{N). We will evaluate 

0{N), 0{\ogN) and 0{N^) terms in log^*^"-'. We will show that these terms are universal 
and the leading order term (the 0{N) term) agrees with the (l,^) ZZ brane disk amplitude. 
We first write down the definition of A^""^ . 

_2-{0,-,0,g„=l,0,-,0} 



^(O-inst.) 

N-1 N-1 



, ^N-lY 



y^^N I dx I n n - ^^f^N-l{K 

Jb<Xi<a 
y{0—inst.) « 

^llS^ / dx {det{x - ^N-ir)t-r' e-^^(^\ (5.1) 

*' Xn. 



where 



Jb<\i<a j_2^ 

„ N-1 N-1 

{det{x - <l>N-ir)t-r'-^ = -^^Z^ / n n - >^^f^N-l{Xl, • • • , XN-lf 



Io-^nst. d<t>N-l deijx - 0^_i)^ e-^^^(^^-^) 

/o-.n,.. d<PN-i e-^try(,._,) 
Here (pN-i is an (A^ — 1) x (A^ — 1) Hermitian matrix. 



(5.2) 
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Keeping (13 .Gj) in mind, we can calculate the factor Z^^_^^ 
()5.1|) as follows. 



(7(0— mst.) 



ZiO-inst.) NhN-1 \^^\NJJ NhN 



(5.3) 



The (9(1/A^) correction comes from the instanton contributions. The calculation of /ijv in 
appendix E of Ref.y^ holds for our case. The result is 



1 



(5.4) 



We need (det(x — 0Ar_i)^)j^_™'^*''' outside the cut, x > a. By setting x = formally, we 
find 



(det(x-0^_i)')!J:r'^ 



/ JJ rfAi Aiv(Ai,---,Ajv)'e- 

J b<\i<a 



(iV-l)!/lo/il---/i7V-2 76<A,<a 
-P/V-l(AAf)^ + 7 Pn-2{.XnY + ■ ■ ■ H 7 -Po('^Af)^ 

PN~i{xf + r7v_iPAr_2(a;)^ H h rAr_irAf_2 ■ ■■riPQ{xf. 



(5.5) 



()5.5|) would hold exactly if there was no limitation to the 0-instanton sector. Actually, there 
is a relative 0{1/N) error coming from the instanton contributions in each term in the last 
line of ()5.5j) . Here we ignored these errors because they turn out to only lead to a relative 
0{1/N) correction in the final result. ()5.5|) is further calculated following Ref. |19j. 



(det(a; - 0,v-i)')S:r'^ 
A;W(x,i; 



g(x, 1) 



exp 



2N I da \ogk^^\x,a) 




1 + 



N 



1 



/ A;W(x,] 

VWJ) J (A;(o)(x,l))^ 



exp 



2A^ [ da \ogk^^\x,a) 
Jo 



1 + 



N 



where 



A;W(x,(t) = i (x-s{a) + ^/ (x — s(cr))^ — 4r(cr) j 



^(x, a) = a/(x - s((t))2 - 4r((T). 



(5.6) 



(5.7) 
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These are 0{N^) quantities. ()5.6|) implies that outside the cut 



V}^]{x) = V{x) -2 [ da \ogk^°\x,a). (5. 
Jo 



This would coincide with ()4.7p outside the cut, so that the structure of the cut in ()5.8p should 
agree with that in W{x). This observation leads to a relation 

(x - - 4r(l) = {x- a){x - b), (5.9) 

from which we obtain 

= s, + 2v^. (5.10) 

From (|5.ip . (|5.3p and (j5.6p . we obtain 

^(") = i- / f ^:^V e-^^^JW (l + o(l)), (5.11) 

where hj\; and Vgj'j-(x) are given in ()5.4j) and ()5.8|) . respectively. From this expression, we see 
that log^'-"^ starts with 0{N) and the saddle point integral in ()5.11|) actually gives 0{N), 
0{\ogN) and 0{N°) terms in log^("). 

First, we assume that V{x) is generic without accidental symmetry. We must treat 
separately the case in which V{x) is even. Substituting ()3.11|) and ()3.22|) into the derivative 
of (j5.8j) and using (j5.1(J|) leads to 



V^l'ix) = V'ix) - [ dre'-'^u { ^ \_ ^ + (9(5^) | . (5.12) 

r^au + TcxX 



As mentioned in section 3, W{x) = —^v}ff i^) starts with a term proportional to 2 in 
the scaling limit, so that we are allowed to simplify (|5.12j) as 

V}^j (x) = iya--2rc ^^'^"^ / dr + 0(e'=), (5.13) 

Y u(r) + a^Vc ^xX 

where we keep only the contribution from the one end r = of the integral region. The 
other terms with integer power in e lower than e'^"^ would be canceled in (j5.12|) . d^V^^j{x) 
starts with the 0{6''~^^) term and N ~ so that Nd^V^^j{x) is finite and we can 
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ignore the 0{6^) correction in ()5.13|) . We also perform the change of the integration variable 
■u(r) = y//iw. Then, using (j3.12p and (|3.14|) . we obtain 

dV^fji^) _i -1 ^.^ f\ 1 ,,,,, 
i4 — =£''^2iia ^fc x^: yjtjfxi 4 dw w-' — — (5.14) 

Substituting ()3.13p into ()5.14|) and using the formula ()B.1|) yields 



d( Vi"""^'''''/ Hi ('" + ;y| ' ^'-if'''). (5-15) 

where Q = a^^Vc ^a;* and Pk-i is the Legendre polynomial of degree k — 1. We have specified 
consistently the lower end of the integral region in such a way that it does not contribute to 
the value of the integral. Furthermore, using ()B.2|) . ()B.3|) and ()2.5|) . we obtain 

^^i = e'^+^ua^n^ — - — fi^-^V2r]k{nn'^(). (5.16) 

2fi; 1 

The above calculation that reduces ()5.14|1 to ()5.1(ij) is essentially same as that in appendix 
B in |12j. The lefthand side of ()5.1(j|l is proportional to the leading term in the scaling limit 
of W{x), which is the universal part of the resolvent, so that we obtain 

Xa = Q~^ =ax:^^ri, Xn = inXa- (5.17) 

That is, X = Xn corresponds to C = ^~^/^^^n and x = a corresponds to C = —Vt^^ii^. By 
integrating ()5.16|) over we finally obtain 

wi;}(x) = {^~'^^^ y "" g^'ii^^-^vuin^i-H). (5.18) 

In order to determine the overall factor K~^z/a5 in (j5.18|) . we need a physical input. We 
adopt the sphere amplitude as the physical input. First, we calculate the sphere contribution 
to the free energy of the matrix model ()3.15|) . By performing the change of the variable 
u{t) = and using (jTT^ . (jSH), dSIIll), (fBl]l and (|R4|l . we obtain 



L — I L2J ^1 
p=0 q=0 

i(K-Va^)2S,(7;V+i (5.19) 



14 



where 

vr 1 



2{2k+l){2k-l){2k-3)' 



(5.20) 



If the sphere amphtude is given by 

^(sphere) ^ ^^^-2^fc+i^ (g ^l) 

where dk is a certain universal constant, K~^z/a5 is fixed as 



= J'^. (5.22) 

Then, NV^'^jr{x) is determined as 

/(O) 

where 



ivv;(°)(x) = D,^?;V^+^^;t(fi/^~^C), (5.23) 



2A; - 1 V Sfc ^ ^ V 2A; - 1 ^ ^ ^ 

We are ready to calculate A^"^ given in ()5.11|1 . 



/i^(A;(o)(x„,l))2 V g(x„,l) 



/i;v (A;W(x„, 1))2 V g(x„, 1) ; * V 

(5.25) 

where we have used ()3.22p . We see from ()5.4|) . ()5.10|) and ()5.17p that in the scaling limit 

= 27ry?;e-^'=^»"'^^^*^*(-^) = 27T,/¥; 



2 



1, 



(A:W(x„,l)) 

/ fcW(x„,l) \^ 1 1 



_ i_ 

Thus, noting VL = a~^rc ^x^:, we finally obtain 



(5.26) 



= I J \ 1 ,|^-!-|e-W^^^*-^(^-), (5.27) 

8 y 7rDfct;;t(^„)^„ + 1 
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where Vk{^), ffe(^n) and are given in ()2.7|) . ()2.10|) . ()2.11|) and ()5.24|) . respectively. It 
follows from this expression that A^^^ is indeed universal i.e. independent of the detailed 
structure in the potential of the matrix model. 

Next, let us consider the case in which V{x) is even. Sn vanishes identically in this case. 
We can see that all the equations (j5.12|) ^ (|5.18|) also hold for this case if we set x^: = l^frl 
and replace a with ^. Namely, 

= («:-Va^)^^^l^^^^;V^^^^;c(^^/i-k~) (5.28) 

with f2 = 4a~^ and Xa = \'^- F^^'p'^'^^^'' given in (j3.15|) does not include s„ so that (|5.19p is 
invariant. However, instead of 1)5.211) we must fix pi'^v'^^^^) as 

p{sphere) ^ 24^7 (5-29) 

since due to the Z2 symmetry there are two critical points each of which contributes equally 
to the free energy. Hence, /t'^z/a^ in (|5.28|) is determined as 

fi;-Va2=2W^. (5.30) 

This indeed reduces fl5.28|) to ()5.23p . ()5.25|) and ()5.26p also hold for this case if we set 
= l^fTc and replace a with |. Noting f2 = 4a~^, we can easily see that A*^") is indeed 
given by ()5.27p . 

Thus, the proof of the universality of A*^") is completed. Note that A^^^ is pure imaginary 
for odd n and real for even n. This reflects the fact that ^ = with n odd corresponds to 
a local maximum and ^ = with n even corresponds to a local minimum. 

As a check on our calculation, let us calculate A^^^ in the k = 2 case, which would coincide 
with /i in The normalization of F^^P'^^'"^) in PI corresponds to ^2 = Using (IT7|l . 
fl2.1()|l . ()2.11|1 and ()5.24j] . we can calculate the quantities that appears in ()5.27j) with k = 2 
and = 1 as 

^2 = -^, ^1 = ^, v,i^,) = -^, v'^i^,) = VQ. (5.31) 
By substituting these quantities into ()5.27|) . we obtain 

^« = \ r/J/x-ie-^sr Vt . (5.32) 

8-3tv^ 
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This indeed coincides with /i in |19j . 

Finally, let us see that the leading order of log^^") indeed agrees with the (1, n) ZZ brane 
disk amplitude, which we denote by Z„. The leading order of log^*^"^ is given by 



-Dug-^^^'^+'^Vk{in). (5.33) 
We can evaluate Zn in Liouville theory by using eqs.(B.4) and (B.6) in [20] as 
, ^,„4-2t (2A;-l)3r(fef) / . 27r V • 2™ ^ 

We can also calculate the sphere amplitude in Liouville theory by integrating twice the two- 
point function of the cosmological constant operators which is given in (2.26) in ^ni- The 
result corresponds to^ 

1 r^2fc-3 . 27r 

dk = ^a/-^ sm— — 7^-4-] ■ (5.35) 

Trsy 2A;- 12A; + 1 2k-l\T{^^)j ^ ' 

It is easy to verify that plugging ()5.35|) into (j5.33p actually yields ()5.34p . 

The fact that the leading order of log^*-"^ is the (1, "n.) ZZ brane disk amplitude is already 
pointed out in • What is new in this section is that we showed that both the leading order 
term (C(l/(7s)) and, in particular, the next to leading order terms ((9(log(?s) and 0{(ff)) in 
log^*-"-* are universal. We also showed that the normalization of the ZZ brane disk amplitude 
is also reproduced precisely by matching the sphere amplitude in the matrix model with that 
in Liouville theory. 



6 2-instanton sectors and the annulus amplitudes be- 
tween the ZZ branes 

In this section, we calculate The estimation of the order in the expansion 

proceeds in the same way as that in the calculation of J^'^^ . We will not dwell on it in this 
section. log^("'"') also starts with 0{N). We will evaluate 0{N), C(logA^) and 0(A^°) 
^The relation between our cosmological constant /i and the cosmological constant /zl in J2| is ^ = 

p/ 2fc-3 \ 

^J■L''^JT^ziZ\■ Note also that for odd k the overall sign of dk obtained from (2.26) in ^] is different from that 

in (|5.35() . It seems possible to attribute this difference to the ambiguity of the sign of the norm in nonunitary 
models. We adjust the overall sign of in (|5.35|) to minus in such a way that it is consistent with the result 
in the matrix model. 
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terms in log . We will see that these terms are universal. We will also show that in 

the n ^ n' case the leading order term in log^i"'" \ which is (9(iV°), reproduces the annulus 
amplitude between the and (1,?^-') ZZ branes, where ^i""'" = A^^'^'^ A'^'^'' A^^'\ First, 

we consider the n ^ n' case. In this case is given by 

2'{0,-,0,g„=l,0,-,0,g„,=l,0,-,0} 



^{0-inst.) 

1 r r r ^'^ 

Zio--I^^(^ -I) dx dy n 

^ J x„ J x„i Jb<\i<a •_i 



x^i Jb<\i<a 



N-2) 



\2 



ry{0—inst.) 
^{O-inst.) 

^^-NV{x)-NV(y)+log{x-y)^ ^ j^g 



N(N-^)^^thn dx rfy(det(x-0;v-2f det(y-0^_2)')!;:r*-^ 



where 



„ N-2 

J b<\i<a ^_-|^ 

(det(x - <PM-2f det(y - <PN-2f)f"''''-^ 



' N~2 
N~2 /N-2 



7(0 



AAr-2)^ 



/o-mst. '^</'A^-2 det(a; - (t)N-2Y det(?/ - 0Af-2)^ e 



N-2 ■ 

2Hptr7, - .^2 „-ivtry(<Aiv-2) 



(6.2) 



Here 0Ar_2 is an (A^ - 2) X (iV - 2) Hermitian matrix. The factor ^JJj""*-)/^(o-mst.) ^j^g 



last line of is calculated as 

y{0—inst.) 

'N-2 ^ 



^(0-inst.) ^ ^2 



Z(0-inst.) JY^^ _ l)hN_^hN-2 N{N - l)h% ' 

(det(x — 0Ar_2)^ det(?/ — (t)N-2Y)%-2^^'^ is evaluated as follows. 
(det(x - (l)N-2f det(?/ - 07v_2)^)!°~"''*''' 



(6.3) 



^tr log(a;-(/)jv-2)^+tr \og(y-<j,N-2f \ (O-inst.) 

I r ' 



I N-2 

(O-m 
'N-2 



exp 



(trlog(a: - ^N-2r)tT"^ + (trlog(i/ - <PN-2r)tT"^ 

+i((triog(x - <PN-2mt-2f + i((triog(y - ^N-2)r)V2r-^ 



+ (tr log(x 

v2\ (O-mst.) , 



(det(x - (f)N~2y)N~T''\^^^iy - 0iV~2)^)!v-2 



X exp 



(trlog(x - 0Ar_2) trlog(?/ - 0Ar_2) ) 



2\(0—inst.) 

)n- 

2\{0—inst.) 
-2 

2\{0-inst.) 



N-2,c 



+ 



(6.4) 



2\ (O—inst.) 



where the subscript 'c' stands for the connected part. We can calculate (det(x— 0Af_2)^))v- 
in a way similar to and ()5.(jp as 



(det(x-</)^„i)')!"-^ 



'Af-2 



2A^ /" rfa logA;(°)(x,a) 



(6.5) 



A similar expression holds for {det{y — 4>N-2f')%-2^^'^ ■ In the leading order of the 
expansion, we can set 

(trlog(x - 07V-2)^trlog(y - (t)N-2f)t~-2T^ = 4(trlog(x - 0)trlog(?/ - 0))^. (6.6) 

The righthand side is calculated in appendix B. From 1)6.11) . ()6.3|) . ()6.4|) . ()6.5|) and ()6.6p . we 
obtain 

^2 . . / wn^ . , X N 2 



A{n,n') _ 



X exp 



/?2 



/ A;W(x,l) 



1 



AW(i/,l) 



V g(x,l) ; (fcW(x,l))4 V ; 



iVV;7;(x) - Wi;j(y) + log(x - yf + 4(tr log(x - 0)tr log(y - 0)), (6.7) 



By noting that in the scaling limit 



(6.8) 



and recalling the calculation of A^'^\ we find 

A{n,n') _ A{n) A(n') A{n,n') 
,yv — "^^c ' 

^(nX) = exp [4(trlog(x„ -0)trlog(x„/ -0))c + log(x„ -a;„/)2 -logrj . (6.9) 

Here ^c"'" is interpreted as the 'connected part' of since A^^^A^'"'"' is the product of 

the 1-instanton contributions. As seen in the previous section, the leading order of log^^"-* 
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is the (1, n) ZZ brane disk amplitude, so that the leading order of log^i"'" ■* is expected to 
be the annulus amplitude between the (1, n) and (1, n') ZZ branes. In the following, we will 
show that this is indeed the case. Using ()B.6|) leads to 

= 4(trlog(x„ - 0)trlog(x„/ - 0))^ + log(x„ - - logr^ 

= log (2xnXn' - (a + b){xn + Xn') + + - 2a/(x„ - a)(x„ - b){xn' - a){xn' -b)^ 
- log (V (.„-a)(.„-6) + V (:^„, -.)(.„, -6) )' 

+ log (^\/Xn - a + \^Xn-b^ + log - a + a/ Xn' - b^ - \og{a - b)^ - log 16 

+ log(x„ -logTc (6.10) 

Recalling that in the scaling limit, 

a = x^{l - Xay/Jj-e), 

Xn = X^{1 + Xa^ny/Ji£), 

we calculate the quantities that appears in ()6.10p . 



(6.11) 



2XnXn' — {a + b){Xn + X„/) + + 6^ — 2^J {xn — a){xn — b){xn' — a){xn' — b) 
= (x,- 6)2(1 + 0(£)), 

a/ {Xn - a){Xn -b) + a/ {Xn' - a) {Xn' - b) = a/x* - 6 {^/Xn - a + ^ Xn' - O) (1 + 0{e^)), 

y/Xn - a + \/xn~b = a/x* - 6 (1 + 0{e^)), 
a/x„/ - a + \/xn' -b = a/x* - 6(1 + O {e^)), 

a-6= (x* -6)(l + 0(e)). (6.12) 
Substituting these into ()6.15|) yields 



log 4"'"') = log 



(x„ - x„0^ {x* - bf 



{^JXn - a + Vx„/ - aY 16rc 
Furthermore, by using ()6.1ip and a relation 

(x, - bf 



16 

which follows from ()5.9j) . we finally obtain 

log^l"'"') = lo£ 



(v/6r+T+ ve^T+T)^ 
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(6.13) 



(6.14) 



(6.15) 



This is universal. The model-dependent quantities such as Vc are indeed canceled. It follows 
that ^("■'"') is also universal because A^"^^ and A^^'"^ are universal. 

Let us see that log^("'"') with n ^ n' is the annulus amplitude between the ZZ branes 
in Liouville theory. First, we rewrite by 



= — sin — = =F'\/ • (6.16) 

2fc-l V 2 ^ ' 

Using this, we express the annulus amplitude between the (1, n) and (1, n') ZZ branes (I2.12j) 
in terms of 

= log- ^-^~^" '^' -. (6.17) 

This indeed agrees with log^i"'"''. 

Next, we consider the case in which n = n'. The same calculation as the n ^ n! case 
leads to 

= ^(^("))2exp [4(trlog(x„ - 0)trlog(x„ - 0)), + log(a;„ - x^f - logrj 

That is, vanishes and log A^"''"'^ diverges. This is consistent with the result in Liouville 

theory where Z„ „ also diverges.^ 

7 Summary and discussion 

In this paper, we analyzed the fc-th multicritical region of the one matrix model that cor- 
responds to the (2, 2k — 1) minimal string theory. We divided the partition function of 
the matrix model in terms of the instanton numbers. We evaluated the ratio of the 1- 
instanton sector to the 0-instanton sector, A^'^\ and the ratio of the 2-instanton sector to 
the 0-instanton sector, We found that log^*^"^ and log ^'^"''^'^ start with 0{N) terms. 

We calculated the 0{N), 0{\ogN) and 0{N^) terms in log^(") and log^("'"'), which cor- 
respond to 0{l/gs), 0{\oggs) and O^g^) terms in the continuum limit, respectively. Our 
results are as follows. {i)0{N), C(logA^) and 0{N^) terms in log^'-"-' and log^*-"'"'^ are 

•^It is already pointed out in 2D^ that the double zero of e'^"" stems from the Vandermonde determinant 
of the matrix model. 
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universal i.e. independent of the detailed structure in the potential of the matrix model. 
(ii)The 0{N) term in log^*^"^ is equal to the {l,n) TIL brane disk amplitude, which is pro- 
portional to (iii)When ^("•"') with n ^ v! vs, expressed as = ^("M''"'^^c"'" \ 
log^i"'" ^ starts with (9(A^°) term. This term reproduces the annulus amplitude between the 
(1, n) and (1, n') ZZ branes in Liouville theory. (iv)The 0{N) term in log^*^"'"^ are given by 
that in log^*^"-* while one of the 0{N^) terms in log^*^"'"^ diverges. This makes log^*^"'"^ 
vanish and is consistent with the result in Liouville theory. 

The above results allow us to assign Figs. 2(a) and 4 to the leading order of log^*-"-* and 
log^i"'"\ respectively. We express ^("'"^ as A^"''^^ = ^{A^"'^)'^A'i'''^\ Then, it is natural to 
assign Fig. 2(b) to the next to leading order (0{\ogN) and 0{N'^)) of log^*^"-' and Fig. 3 to 
the leading order of A^'^\ which is divergent. Namely, the next to leading order of log^*-"-* 
can be interpreted as the annulus stretched within a single (1, n) ZZ brane while the leading 
order of log can be interpreted as the annulus stretched from one (l,n) ZZ brane to 
the other (l,n) ZZ brane. In Liouville theory, these two diagrams are not distinguished 
because the ZZ branes have no intrinsic parameter such as the position except the label n, 
so that the annulus amplitude between the two identical (1, n) ZZ branes in Liouville theory 
is consistently divergent. As is stressed in ^U] in the k = 2 case, a nontrivial thing is that the 
next to leading order of log A^"^^ is a finite and universal quantity, which cannot be evaluated 
at least so far in Liouville theory. A^'^^ has a physical interpretation as the chemical potential 
of the n-th instanton, so that it is suggested that the matrix model possesses the information 
on the nonperturbative effect that Liouville theory cannot predict. Note that as shown in 
[in] this quantity cannot be calculated through the loop equation (the string field theory), 
either. 

Our results imply that vanishes if g„ > 2 at least for a certain n so that the 

expansion of the partition function is terminated with 2^~^ terms as 

The maximum of the total instanton number is q = k — 1. Actually, A^^'^''"'^^ is nonva- 
nishing. Our results suggest that each (j7.ip is universal and its 'connected 

part' systematically reproduces the amplitudes among the ZZ branes in Liouville theory. It 
is interesting to calculate multi-point amplitudes among the ZZ branes in Liouville theory 
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if possible and compare them with ^{'Jir-^'jfc-i} The reason why with 

g„ > 2 at least for a certain n vanishes is that one cannot place more than one eigenvalue at 
the same point because of the repulsive force coming from the Vandermonde determinant. 
However, taking into account log(x„ — Xn')"^ stemming from the Vandermonde determinant 
from the beginning, one can obtain a 'classical solution' that realizes the situation in which 
those eigenvalues are separated a little from each other around a local extreme of the poten- 
tial. Perhaps one can construct an expansion of the partition function based on this classical 
solution, which differs from the expansion in this paper. This expansion would lead to the 
singularity destroying deformation argued in ^Dj, which changes the singularity to a cut, 
and may yield (9((e^")™') (m > 2) corrections to ()7.1|) . Note that in the k = 2 case one can 
reproduce from (jT.lj) with these possible corrections the prediction of the string equation 
that the deviation from the perturbative solution for the free energy F = log Z behaves as 
~ e^^ at leading order of e^^, because at this order the logarithm of Z in ()7.1|) with the 
possible corrections takes the form F = pi^-'^"-'^*-) -f 

Finally, we make a comment. As is pointed in [171120]! the sign of NV^^j-{^n) is (— 1)""*"^ 
so that ^'•"^ with n even behaves as ~ e'^^ , which is catastrophic. The energy of n-th 
instanton with n even which corresponds to the local minimum is below the Fermi level. 
Therefore, the perturbative vacuum is unstable due to the eigenvalues tunneling to these 
local minima. This is due to the nonunitary nature of the model. Note that the fact that 
in the (2, 3) case n takes only 1 is consistent with the unitarity of the (2, 3) model. Thus, 
the expansion of the partition function of the one matrix model in the instanton numbers 
should be understood formal one in this sense. 

It is important to generalize our analysis to the two matrix model, which can represent 
the unitary noncritical string theories. 



(a) (b) 



Figure 2: log^^"^: (a)the leading order (b)the next to leading order 
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Figure 3: The leading order of log^c 
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Appendix A: The (2, 5) model 

In this appendix, as an example, we analyze in detail the one matrix model in the 3rd critical 
region {k = 3), which corresponds to the (2,5) minimal string theory. We adopt an even 
potential of the one matrix model, 

V{x) = ^x^ - ^g^x' - ^(76x6. (A.l) 

The first equation in ()3.8j) gives in the large- iV limit 

a = r{a) - 3g,r{af - Wg^riaf, (A.2) 

while the second one is trivial. Correspondingly, s(cr) vanishes identically. The scaling limit 
fj3.11|) is given in this case by 

96 = -^(l-3/5/i£^), 
24 



r{a) = 3 ^1 — -aeu^r] 

(7=1 — e^UT, 

— = e2Kgs. (A.3) 



Plugging ()A.3|1 into ()A.2j) gives 



3/3q; a 



3 



/iM + — M^. (A.4) 



2i/ " 8z/ 

Comparing this equation with the string equation ()3.12|) . we find 



8u V8' 2u 

which are equivalent to 



1 



(A.5) 



a = 2^3/3^, u=^3V3p^. (A.6) 

W{x) takes the form 



W{x) = -ge{x^ - xl){x^ - xl)Vx^-a'^. (A.7) 



The condition that R{x) ~ ^ when a; — > cxd is equivalent to 

94. + 96 + xl + xlj = 0, 

1 I ^ i ^ 2/ 2 , 2\ 2 2 

1 - 5-6 I gCt - -a (Xi + X2) - 

-^76 - + xl) + ^a'xlxl^ = 1. (A.8) 

In the scaling limit ()A.3jl . a, xi and X2 are determined by these equations as 



X2 = X. |l + ^i±^v^el, (A.9) 
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where = 2-\/3. If we put x = + eQ, we obtain 

5 



\/2 5-5 



540 



where 



= — , 

Xa 

1-V5 2tt 

ti = = — cos — . 

^ 4 5 ■ 

1 + ^5 Air 

to = = — cos . 

^ 4 5 



By integrating W{x) over x, we obtain 



NV^ff{x) = -2N dx'W{x') = ^g~^K-^xlxl^~'^iiv^{n^i-H) 

4 ■ 3 4 7 7 1 ~ 

It is easy to verify that the last hne in ()A.12|) can also be obtained by plugging k 
(EH) into (jSI2Hl)- 

Appendix B: Useful formulae 

In this appendix, we gather some formulae, which we use in sections 5 and 6. 
The following formulae for the Legendre polynomials are used in section 5. 

\-\ 

^ ^ ^ (-ir(2n-2p)! 
" (n-2p)!p!(n-p)! 

Pn{-x) = (-irP„(x). 

1 Tn{x) + Tn+l{x) 



/X 
dt {X - t)-^Pn{t) 



n + \ y/x+l 
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7 mn / ^ m(m - 1) ■ ■ ■ (m - n + 2) 

dxx"^PJx) = - ^ ^- -. (B.4) 

[m + n + l)[m + n — 1) ■ ■ ■ [m — n + 3) 

The cylinder contribution to the two macroscopic loop correlators in the one matrix 
model with a general potential is obtained in [23 • The result is 



1 \^ f 1 ^\ 1 / xy — ^{a + b){x + y) + ab 



^'[x-<p)^'[y-<p)/^ 2{x-yf {^^ - a){x - h){y - a){y - h) ^^'^^ 
By integrating this over x and we obtain 

(tr log(x - 0) tr \og{y - 4>))c 

dx' I dy' / tr ( — ^ — ) tr 



X J y 



- log 

2 ^ 



x' — (j) J \y' — 4>, 

2xy - {a + b){x + y) + a"^ + b'^ - 2a/(x - a){x - b){y - a){y - b) 
- log (y{x-a){x-b) + a/(z/ - a){y - b)^ 
+ log {^x — a + \/x — bj + log ^Vz/ ~ + \/y — bj — \og{a — b) — log 2. (B.6) 
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